Three-body problem in Fermi gases with short-range interparticle interaction 
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We discuss 3-body processes in ultracold two-component Fermi gases with short-range intercompo- 
nent interaction characterized by a large and positive scattering length a. It is found that in most 
cases the probability of 3-body recombination is a universal function of the mass ratio and a, and 
is independent of short-range physics. We also calculate the scattering length corresponding to the 
atom-dimer interaction. 
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In the search for fermionic superfluidity, ultracold 
Fermi gases with strong interactions and high densities 
are produced routinely. In current experiments with 
[1] and ^Li [2-5], the scattering length a, correspond- 
ing to the interaction between different fermionic species, 
has been modified using a powerful tool of Feshbach res- 
onances. Reaching high densities is limited by 3-body 
recombination - the process in which two atoms form a 
bound state and a third one carries away the binding en- 
ergy e. In the case of a two-component Fermi gas the 
3-body process requires at least two identical fermions 
to approach each other to a distance of the order of the 
size of the final bound state. Therefore, in the ultra- 
cold limit the recombination probability acquires an ad- 
ditional small factor K/e, where K is the relative kinetic 
energy of identical fermions (cf. [6]). One may speculate 
that limitations on achieving high densities are not as 
severe as in Bose gases, where the recombination rate is 
independent of the kinetic energy of particles. 

Theoretical studies of the 3-body problem have re- 
vealed the existence of two universality classes for the 
case where a greatly exceeds the characteristic radius 
of interatomic interaction. Re [7,8]. For the first class, 
short-range physics is not important and the 3-body 
problem can be described in terms of 2-body scatter- 
ing lengths and masses of particles. One can then use 
the zero-range approximation for the interatomic poten- 
tial, which has been successfully demonstrated, e.g. in 
the calculations of neutron-deuteron scattering with the 
total spin S — 3/2 [9]. For the second class, where short- 
range behavior is important, the description of the 3- 
body problem requires at least one parameter coming 
from short-range physics [10]. For two identical fermions 
(t) interacting with a third particle (|), the presence of 
two universality classes has been shown by Efimov [7,8]. 
For a large a, the pair interaction between distinguish- 
able particles leads to the appearance of an attractive 
interaction in the 3-body system [11]. If the mass 
ratio m|/TO| is smaller than approximately 13, then this 
attraction is not sufficient to overcome the centrifugal 
barrier l{l + l)/r^ = 1 due to the symmetry). In 



this case the probability of all three particles being in 
the volume vanishes and short-range physics drops 
out of the problem. However, the effective attraction 
increases with the mass ratio and for m| /to| ^ 13 we 
have a well-known phenomenon of the fall of a particle 
into the center in an attractive potential [12]. The 
shape of the wavefunction at distances of the order of i?e 
then significantly influences the large-scale behavior and 
short-range parameters of the interaction potential are 
required to describe the 3-body system. 

In this paper we discuss a 3-body system (tti) con- 
taining two identical fermions and belonging to the flrst 
universality class. In the ultra-cold limit we find uni- 
versal functions of a and m| /m| for the probability of 3- 
body recombination and for the amplitude of atom-dimer 
scattering. As expected, the recombination to compara- 
tively deep bound states is much slower than in the Bose 
case. However, despite the suppression factor of K/e, 
the recombination to a weakly bound s-level (a > and 
a ^ Re) for realistic parameters of a two-component 
Fermi gas can be as important as in a Bose gas with the 
same density and scattering length. 

In the center of mass reference frame the state of the 
3-body system (t t i) with total energy E is described 
by the Schrodinger equation 
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where hy/y^rnj is the distance between identical t" 
fermions, hx/^/2fl is the distance between their cen- 
ter of mass and the |-particle, /x = 2m^mi/{2m^ + 
mj^) is the corresponding reduced mass, and 6 ~ 
arctan \/l + 2r7i|/m|. The vector X — {x, y} describes 
the rescaled 6D configuration space for the 3-body prob- 
lem. The potential of interaction between distinguish- 
able particles is V and the interaction between identical 
fermions is omitted. 

We assume that the range Re of the interatomic po- 
tential V is much smaller than the (positive) scatter- 
ing length a, which in the rescaled coordinates equals 
l/\/e. Therefore, in the 6D space the rhs of Eq.(l) van- 
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ishes everywhere except in the close vicinity to 3D sub- 
spaces S+ and S-, which we define parametrically as 
X±{^) = {rcos^, ±r sin^}, where the 3D vector r is a 
parameter. Near the subspace 5+ (S*-) the J,-particle is 
close to one of the t-fermions and is separated by the 
distance r from the other. The distance to the first t- 
fermion wc denote by r^, which is then a parameter in 
a 3D subspace orthogonal to S+ {S-). The subspace 
S- differs from S+ by the permutation of the f-fermions 
(y — > — y). We thus face a typical boundary- value prob- 
lem for the Poisson equation [— V^f — iS] \l/ = 0, with a 
boundary condition given by 



*~ 1- l/(v^|r^|) + 0(rj.). 



(2) 



Introducing functions /+(r) and /-(r) we can write 
a general solution of Eq.(l), which is valid everywhere 
except the narrow vicinity of S+ and S- : 

^ = '^o{X) + J2 [Ge {\X - X±(r)|) f±{r)cfT. (3) 



Here '^q is a solution of the Poisson equation and has 
no singularities at S+ or S-. We restrict 4*0 to be finite 
everywhere and, hence, for £ < we have = 0. The 
Green's bmction is a solution of Eq.(l) with the rhs 
5^{X) = (5^(x)(5^(y) and is given by 



~EK2(^/^\X\) 
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E <0 
E>0 
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dimer binding energy e. The total energy of the system 
is then negative E !v —e and, consequently, "ifo = 0. As 
the size of the dimer is much smaller than the inverse 
relative momentum, the s-wave contribution dominates 
and Eq.(4) becomes 



L'=7(r) = 0. 



(6) 
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Here is an integral operator which is obtained from 
L_£ by integrating over angles. 

Prom Eq.(3) we find that far from the origin (r ^ 
\l\fe) and in the region of S± the wavefunction 
±/(r) exp(— •\/£r_L)/47rr_L. Therefore, at these distances 
the function /(r) describes the atom-molecule relative 
motion and behaves as 1 — P/r. The atom-molecule scat- 
tering length is then given by am = af3y/s sm20. The 
ratio am /a is plotted in Fig.l as a function of m^/mi. 
In the limit of m| / 3> 1 one can use the Born- 
Oppenheimer approximation. In this case the heavy f- 
fermions move slowly in a field produced by the exchange 
of the fast light J,-particle. The adiabatic behavior as- 
sumes decomposition of the wavefunction into two parts. 
The first part describes s-wave scattering of the heavy 
fermions. The motion of the light particle is, therefore, 
described by a wavefunction antisymmetric with respect 
to their permutation. At large distances (y » the 
effective interaction has the form of a repulsive Yukawa 
potential U{y) = cos0{mi/mi)2^/£exp{— 2 cos 9 y-/e)/y. 
The corresponding scattering length is plotted in Fig.l 
as a dashed line. We estimate am/ a ^ ln(m^/mj^). 



Here K2 is an exponentially decaying Bessel function and 
H2 is a Hankel function representing the outgoing wave. 

In order to satisfy the symmetry with respect to per- 
mutation of identical fermions we set /+(r) = — /_(r) = 
/(r) and require 5'o(x,y) = — 4'o(x, — y). Eqs. (2) and 
(3) then give the integral equation for the function /(r) 



4nMX+{r)), (4) 



where the hermitian operator Le is given by 
W(r) = 47r y [Gb (|r - r'l) [/(r) - /(r')] 

+ Ge (\/r2-hr'2-2rr'cos26l) /(r') 



d'r'. (5) 



This operator conserves angular momentum, and we can 
expand the solution of Eq. (4) in spherical harmonics and 
deal only with a set of uncoupled integral equations for 
functions of a single variable r. Eq.(4) is a particular 
case of the Faddeev equations [13] and in the momentum 
representation it was first obtained in [9] for m-^ = rn^ . 

The method presented above allows us to solve the 
problem of atom-molecule scattering for an arbitrary 
mass ratio. Let us assume that the relative kinetic en- 
ergy of a ti-dimer and |-fermion is much less than the 
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FIG. 1. The ratio am/ a versus m-f/mj. calculated from 
Eq.(6) (solid line) and from the Born-Oppenheimer approxi- 
mation (dashed line). 

It can be shown analytically that /(r) ~ r'^ near the 
origin. The exponent 7 increases monotonically from 1 
to 3, with m^/mi growing from to 00 indicating an 
effective repulsion in the system [7] . In the case of m-f = 
mi we obtain am/a ~ 1-2 (cf. [9]). 

We now turn to the problem of 3-body recombina- 
tion just above the threshold < <C e. This pro- 
cess (T + T + i — > T + Ti) takes place at interpar- 
ticle distances comparable with the size of a tJ,-dimer 
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l/\/e). Therefore at these distances the shapes of the 
functions \E' and / are independent of E. The normal- 
ization coefficient, however, depends on the large scale 
behavior and, therefore, on energy E. We can illustrate 
this by writing the wavefunction for free motion in the 
form 4'o(x, y) = sin(k-y) exp(iq-x), where the momenta 
k and q satisfy the equation E = + . At distances 
comparable with the function \l/o reduces to k • y, 

giving rise to the threshold law for the probability of 
3-body recombination [6]. Considering the limit E ^ 
we keep the term 47rk • r sin 9 in the rhs of Eq. (4) . 

We now search for the solution of Eq.(4) in the form 
/(r) = /(r)k • r/kr, which corresponds to the angular 
momentum 1 = 1. Eq.(4) then transforms to 



(L[j=i - Vff) fir) = Ankrsine. 



(7) 



An explicit form of the operator Z[f ^ can be obtained 
from Eq.(5) by integrating out the angular part of /(r). 

The operator L'lp^ has remarkable properties which al- 
low us to solve Eq.(7) analytically. Firstly, for /(r) = r'', 
where v is in the region —4 < Re (ly) < 2, the action of the 
operator reduces to Ljp^r'^ = X{i')r'^~^ , and the function 
X{u) is given by 



i^cos26>cos [{ly + l){2e - 7r/2)] + sin [t/(26l - 7r/2)] 
l)cos2 26'sin2^sin(7riy/2) ' 

In the specified interval of v the function A has two roots: 
uq and —1/0—2. With m-f / mi increasing from to approx- 
imately 13.6 the root decreases from 1 to -1 approach- 
ing the second root. At larger mass ratios the roots are 
complex conjugates with the real part equal to -1. Zeros 
of A(z/) are important to understand the behavior of /(r) 
at short distances. 

Secondly, any eigenfunction x{f) of the operator Lq^^ 
corresponding to the eigenvalue 1 can be rescaled to gen- 
erate an eigenfunction corresponding to a new eigenvalue 
p (for any p > 0) 



(8) 



As A has two roots, the spectrum of L\p^ is two-fold 
degenerate, and for p = 1 there is only two linearly in- 
dependent eigenfunctions which behave as xi ~ r'^" and 
X2 ~ T-^'^o-a small r. The complex exponents i^o and 
Vq at larger mass ratios imply that the function xi = X2 
rapidly oscillates as r ^ and the Efimov effect takes 
place. In both cases we choose two real eigenfunctions 
which at large distances (r ^ 1) have the asymptotes 
Xi,2('') ~ 2sin(r -|- di,2)/r. The phases and 62 can be 
determined from a numerical calculation of Eq.(8). 

The case p = ^/e corresponds to the collision of a 
molecule with an atom just below the 3-body recombina- 
tion threshold {E = —0) where the function ^0 is zero. 



Strictly speaking this process is described by the linear 
combination fmoi{r) = Axi{y/sr) + Bx2{V^r). The de- 
termination of the coefficients A and B involves short- 
range physics, i.e. the knowledge of the 3-body wavefunc- 
tion "if at distances of the order of Re, and is beyond the 
scope of this paper. However, the matching procedure 
implies that at these distances both terms in fmoiii") are 
of the same order of magnitude, so B/A ~ {^/eRe)^'^°~^'^. 
In the case of mi /mi < 13.6, the exponent 2fo -|- 2 
is real and positive, and at distances r ^ Re one has 
fmoii'r) ~ Axi{\/sr). Hence, short-range 3-body param- 
eters are unnecessary. 

Fortunately, the discussed case of m| /m^ < 13.6 covers 
almost all practical situations. Using the function Xi(r) 
and the scaling property (8) we can construct a com- 
plete orthonormal set of functions /p(r) = PXi{P''")- Com- 
pleteness follows from the equality rfp{r) = pfr{;p) and 
from the normalization condition fp{r)fp'{r)r^dr = 
2tiS{p — p'). The operator Lq^^ — \/e in Eq.(7) can now 
be inverted and after some manipulation we obtain 

47rfcsin6» ( ^ A(0)A(1) 
/(r) = + A(l) + ^ ' ^ ' 



A(0)A(1) /A(-l) 
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The integral in Eq.(9) is taken as a principal value. The 
solution / is not singular at short distances. The terms 
containing r^^, r° and r vanish as expected except in 
situations where A(l), A(0) or A(— 1) are equal to 0. At 
distances r <C l/-\/e we have /(r) ~ and the short- 
range physics does not come into play. For large distances 
the asymptotic behavior xiix) k, 2sin(r -|- (5i)/r gives 
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^^^")"^"=2.cos(Vir + ^r) + 0^ ' 
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Obviously, this describes the molecule-atom channel re- 
sponsible for 3-body recombination. An outgoing wave 
of the form exp(iY^r)/r describing the dimer and atom 
flying apart, is obtained by adding a general solution of 
the homogeneous form of Eq.(7), which is proportional 
to fmoi{r), to /(r). 

Eqs. (9), (10) and (3) provide the 3-body wavefunc- 
tion ^{X) at distances < \X\ < 1/VE, where 
the inelastic dimer-atom channel is well separated from 
the elastic one, and give the amplitude of recombina- 
tion. Then for the number of recombination events 
(T + T + i — > T + Ti) per unit time and unit volume in 
a gas we obtain 



^lll = ai{ei/e)nin 



V 



(11) 



where and are the densities of the fermionic com- 
ponents. The rate constant a^{€^/e) is proportional to 
the average kinetic energy of t-particles , which equals 
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3T/2 in a nondegcncrate gas, and 3Tp/5 in a deeply de- 
generate gas. The coefficient is given by 

= ^-£2A2(l)A^(0)A(-l)sin2^tan3^. (12) 

Fig. (2) shows the dependence of the dimensionless quan- 
tity a-^m^e^ /h^ on m-^/mi. In the case of m-f = mj^ = m 
the quantity a-f « 148/ia^/m. 
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FIG. 2. The parameter a-^m^e^ /fC' versus m|/m|. 



The absence of 3-body recombination at the points 

A(l) = {mi^/m^ = 0), A(0) = (mt/mx w 8.62) and 
A(— 1) = {mi /mi « 13.6) is a purely quantum phe- 
nomenon. Prom Eqs. (9) and (10) we see that for these 
mass ratios the free atom channel is decoupled from the 
molecular recombination channel and the interparticle in- 
teraction leads to elastic scattering only. 

The process (i + i + t — > i + it) is described in a 
similar way. The results are given by Eqs. (11) and (12), 
where one should interchange the subscripts t and J.. In a 
mixture of two hypcrfinc components of the same isotope, 
with equal densities = ni = n/2 and ef = ej = e, the 
total loss rate of particles is 

-n/n = Ln^ « 111 {na^f e/h. (13) 

This formula is valid provided the inequalities e ^ 
fi' jmc? <^ fi^ /m,Rl are satisfied. In current experiments 
with and ^Li one has e ~ 15/iK, n ~ lO^'^cm"'^, and 
the scattering length can be tuned using Feshbach res- 
onances. For a « 280A we calculate the rate constant 
L - 10-25cm'^/s and the loss rate (13) is ~ 10 s"^ 

In a trap with a barrier etr higher than the kinetic 
energy of dimers and fast atoms produced in the re- 
combination process (~ Ti? /ma^), these particles remain 
trapped and Eq.(13) overestimates the losses. Actually, 
the peak loss rate is reached at a finite scattering length 
a ^ y^mitr//i, not for a ^ oo. This is consistent with re- 
cent observations of particle losses in trapped ^Li, where 
a was tuned by a Feshbach resonance [2-5]. 

We emphasize the difference between two-component 
Fermi gases and Bose gases with respect to inelastic pro- 



cesses. Energy dependence of the recombination in Fermi 
gases leads to decreasing e since particles with compar- 
atively high kinetic energies recombine. However, the 
degeneracy parameter decreases due to a faster particle 
loss. Eq.(13) gives an estimate of the 3-body loss rate 
even if a ~ iZg (i-e. deep bound state). In Fermi gases 
the dependence (13) leads to a substantial suppres- 
sion of the recombination probability compared to the 
dependence in Bose gases [14]. Further, in contrast 
to Bose gases, cold collisions between dimers and atoms 
have fewer inelastic channels. Indeed, there are no weakly 
bound Efimov tti-states if m^/mi < 13.6, and the effec- 
tive repulsion between an atom and a dimer at low ener- 
gies suppresses the relaxation to deeper molecular states. 
In the case of a deep trap and large a, where products 
of 3-body recombination stay trapped, one expects a fast 
creation of weakly bound molecules which can be further 
cooled. The lifetime of the gas of these Bose dimers is 
sufficient to observe, for example, a BEC. 
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